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We introduce a semi-classical limit for many-body localization in the absence of global symmetries. 
Microscopically, this limit is realized by disordered Floquet circuits composed of Clifford gates. In 
d = 1, the resulting dynamics are always many-body localized with a complete set of strictly local 
integrals of motion. In d > 2, the system realizes both localized and delocalized phases separated 
by a continuous transition in which ergodic puddles percolate. We argue that the phases are stable 
to deformations away from the semi-classical limit and estimate the resulting phase boundary. The 
Clifford circuit model is a distinct tractable limit from that of free fermions and suggests bounds on 
the critical exponents for the generic transition. 
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The central assumption of statistical mechanics is that 
interactions between particles establish local equilibrium. 
It is now believed that this assumption of ergodicity 
breaks down in isolated quantum systems in the presence 
of sufficient quenched disorder [IHIl]. The fundamental 
feature of such many-body localized (MBL) systems is the 
slowness of the spread of entanglement [mun] due to the 
presence of local integrals of motion (lOM) p [T2ll2Tll^ . 
These local observables remember initial conditions for¬ 
ever and thereby block both equilibration and transport. 
At lower disorder strengths, the system delocalizes and 
ergodicity is restored. There is now much theoretical 
and numerical evidence for this picture in one dimension 
[IHIl [HHH] and the experimental search in quantum 
optical and atomic systems is underway [251E7H29]. 

The nature of the two phases and the dynamical tran¬ 
sition between them is poorly understood in dimensions 
greater than one. In this article, we present a semi- 
classical theory of the delocalization transition in the ab¬ 
sence of global symmetries in arbitrary spatial dimension. 
In the semi-classical limit, the transition maps precisely 
onto classical site percolation. Each lattice site either 
blocks the flow of quantum information or permits it to 
spread. Clusters of unblocked sites form ‘ergodic pud¬ 
dles’, as such regions locally thermalize. Blockage sites, 
on the other hand, host local integrals of motion. Thus, 
at large density of blockage sites the system is lo¬ 
calized. As Ply decreases, the ergodic puddles grow and 
eventually percolate (p^) leading to delocalization. 

The sharp distinction between blocked and unblocked 
sites is a consequence of the lack of quantum tunneling 
in our semi-classical limit. More precisely, the localiza¬ 
tion length of the integrals of motion on blocked sites 
is identically zero. The delocalization transition is in¬ 
stead driven by the divergence of the percolation length 

characterizing the size of ergodic puddles. Perturbing 
away from the semi-classical limit, both length scales be- 
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FIG. I: Schematic dynamical phase diagram without conser¬ 
vation laws. The semi-classical/Clifford limit corresponds to 
the horizontal axis (^g = 0). The blue line covers the delocal¬ 
ized non-ergodic phase ending at the delocalization transition 
at pb = Pb in the percolation class. In d = I, pg = 0. Car¬ 
toon: red regions show ergodic puddles with pink regions of 
influence due to tunnelling. 


come important as quantum information may now tunnel 
through blocked regions. Fig.[^shows the resulting phase 
diagram and the generic delocalization transition. 

The semi-classical model has several desirable features. 
First, the absence of symmetries distills the underlying 
entanglement transition from the complicated epiphe- 
nomena associated with globally conserved currents. Sec¬ 
ond, it is computationally tractable m but the phase 
diagram does not coincide with that of the usual free 
fermion Anderson models [31]. Finally, although the 
phases and transition are semi-classical, they provide a 
new perturbative starting point for studying the generic 
transition. 

The model — The breakdown of ergodicity is simplest 
in the absence of global symmetries. The equilibrium dis¬ 
tribution is unconstrained and there are no slow modes 
associated with conserved currents. We therefore study 
discrete unitary dynamics generated by the periodic ap¬ 
plication of a unitary to a system of N qubits — that is, 
the state after t time steps is \'ip{t)) = t/^|'0(O)). Such dy- 
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namics need not conserve anything, including energy, and 
are known as Floquet dynamics; see [32H34] for recent 
studies of MBL in Floquet systems with time-varying 
Hamiltonians. 

We implement as a multi-layer quantum circuit. 
Each layer is composed of non-overlapping nearest neigh¬ 
bor two qubit gates on an underlying spatial lattice, while 
adjacent layers are stacked in a brick-layer configuration. 
An example in d = 1 is shown at the bottom of Fig. 

In the Heisenberg picture, the support of bounded oper¬ 
ators may only grow in discrete steps under the evolution 
of I/, which provides a technical advantage over generic 
time-varying Hamiltonians. 

The quenched disorder in the model comes from the 
choice of the particular two qubit gate at each brick. Re¬ 
stricting this choice to members of the Clifford group de¬ 
fines our semi-classical limit; we will see below = 0 for 
such circuits. The Clifford group on N qubits consists of 
those unitaries which preserve the group {/, X, T, 
of Pauli operators under conjugation. This is clearly a 
special subgroup of U{2^); generic unitaries map Pauli 
operators into superposition. Moreover, the Gottesman- 
Knill theorem [35] guarantees that the time evolution of 
Pauli operators by such Clifford circuits is classically ef¬ 
ficiently simulable. This provides both analytic and nu¬ 
meric tractability. For reviews of the many properties of 
Clifford circuits, see mmm- 

Global quenches from product states lead to ther- 
malization in ergodic circuits. That is, the reduced 
density matrix of any finite subsystem A approaches 
the completely mixed state, limt^oo limAr^oo 1. 

Equivalently, the entanglement entropy of A, SA{t) = 
—Trp^(t) log 2 (pA(^)) is maximal at late times. This pro¬ 
vides a simple test for ergodicity. A circuit with a local 
integral of motion O fails this test as the local density ma¬ 
trices cannot forget the initial value of O even as t ^ oo. 
We remind the reader that a local integral of motion O is 
a local operator conserved by the evolution: [I/, O] = 0. 

One dimension — One-dimensional systems are espe¬ 
cially prone to localization in the presence of disorder. 
In the Clifford limit, there are local motifs of gates which 
entirely block the spread of operators in the Heisenberg 
picture. A finite density of such motifs thus chops the 
chain into puddles from which information about local 
initial conditions cannot escape, preventing thermaliza- 
tion. 

To illustrate this, consider a restricted ensemble: each 
gate in U is chosen to be either B or R (blue and red ovals 
in Eig.|^ independently with probability p and 1—p. The 
Clifford gates B and R are defined (up to inconsequen¬ 
tial signs) by their action on the Pauli generators under 
conjugation [54]: gate R: ZI ZZ, XI ZY, IZ ^ 

XX, IX YX, gate B: ZI ZI, XI XZ, IZ 

IZ, IX ^ ZX. 

Eigure [^ shows three representative Heisenberg evo¬ 
lutions of Pauli operators under the action of a circuit 
drawn from this gate set. There are several distinct be¬ 
haviors. Within the domain of B gates on the left, lo- 



FIG. 2: Space time evolution of Pauli operators Za{t), Xb{t) 
and Zc{t) under periodic evolution by the d = 1 brick-layer 
Clifford circuit U shown below. At each time t, the sites are 
colored red, blue, or green according to the presence of X, 
Y or Z respectively. Dark (light) qubit lines are blockages 
(non-blockages). Blue (red) ovals denote gate B (R). 


cal operators do not spread. The Z operator beginning 
at site a is conserved by the evolution, representing a 
strictly local integral of motion. The conjugate X oper¬ 
ator at site a behaves just like that shown at site b: its 
spreads to the neighboring sites and then oscillates back. 
Those sites touched exclusively by B gates are examples 
of bloekages. A blockage site hosts an lOM which remains 
on-site after each layer within the circuit U. Elementary 
arguments show that such sites prevent the evolution of 
operators with support strictly on the left (right) from 
extending past to the right (left). 

The Z operator beginning at site c spreads uniformly 
within a light-cone inside the region of R gates on the 
right. This expansion is blocked by the blockage on the 
left, but passes around the isolated B gate on the right. 
Any other local operator beginning at site c exhibits a 
similar expansion behavior (not plotted). Intuitively, 
when local Pauli operators evolve into long strings, we 
expect the reduced density matrix to forget local ini¬ 
tial conditions. Thus, we conjecture that clusters of 
non-blockage sites act as fully ergodic puddles, wherein 
Pa(^ ^ oo) oc IL for subsystems A less than half the size 
of the puddle. We have checked this conjecture numer¬ 
ically for initial product states in which each spin inde¬ 
pendently points along the x,y or z direction, but we do 
not have a general proof, d We emphasize that blockage 
sites are not trivial cuts in the chain, as X and Y are 
not additional integrals of motion. Rather, the support 
of X{t) (Y (t)) on a blockage site grows until other block¬ 
age sites are encountered, as in the oscillations shown by 
X(t) at site b. 

Putting the pieces together, the density of blockage 
sites is pb = p^ for a /c-layer circuit drawn from the re- 
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stricted ensemble. At any > 0, these blockages prevent 
the percolation of quantum information across the chain. 
This Clifford localized phase has an extensive algebra of 
local lOMs of two kinds: the single site lOMs on the 
blockage sites and larger scale lOMs associated with the 
ergodic puddles (of typical size ^e)- Both types of lOMs 
are strictly local {^q = 0) as they have bounded support 
on the lattice. At small the typical puddle size 
diverges as l/p^ and the system is ergodic at = 0. 



FIG. 3: (left) Growth of entanglement entropy of subsystem 
A after global quench from z-product state in square lattice 
Glifford circuit, with A = 40 x 40 and Na = 24 x 24. (right) 
Subsystem size Na dependence of disorder averaged saturated 
entropy at long times [5'a(oo)]. Dashed lines indicate volume 
(slope 1, green), fractal (slope = 91/96, black) and area 
(slope 1 — l/d=l/2, red) power laws. 


entanglement entropy after a global quench (see Fig|^. 
The strict locality {^q = 0) of the lOMs in the MBL phase 
implies that there is no slow growth of entanglement m 
|2Q] . but rather, the entanglement entropy 5 'a(^) saturates 
to an area law [39] as seen at p^ = 0.9. In the delocalized 
phase at p^ < p^, the saturated entropy satisfies a volume 
law (p 5 = 0,0.2). At p 5 > 0, however, the finite density 
of blockages suppresses the coefficient relative to the fully 
ergodic case. Thus, the semi-classical delocalized phase 
is not ergodic. 

At the critical point, Pb = pI ^ 0.4, the intersection of 
the percolating cluster with subregion A contains Np ^ 

jy^f/d where df is the fractal dimension. It follows 

that Sa{oc) . This anomalous growth law at 

the critical point is rather unusual. It is only allowed 
because the delocalized phase is non-ergodic, as predicted 
by Ref. [40j [55] . 



Higher dimensions — The generalization to d > 1 is 
immediate. Sites may host blockages, which prevent the 
spread of operators under Heisenberg evolution. If the 
ergodic puddles of unblocked sites percolate, the circuit 
delocalizes; if not, it is MBL. As usual for site percola¬ 
tion, the percolation threshold p^ is non-zero in d > 2 so 
that both phases arise in the semi-classical limit with a 
transition in the percolation universality class. The as¬ 
sociated critical exponents are well known [37l [38] . For 
example, the characteristic size of the ergodic puddles 
grows as ^ \Pb — Pb\~^ Pb Pb from above with 
z/ = 4/3 in d = 2, u ^ 0.88 in d = 3 and u = 1/2 for 
d > 6, the upper critical dimension. The MBL phase 
shrinks in higher dimensions as there are more paths for 
operators to evade blockage sites. On the hyper cubic 
lattice, the metallic clusters stop percolating at pg ^ 0.4 
in d = 2, p^ « 0.7 in d = 3 and p^ « 1 — l/2d to leading 
order in large d. 

We note that the site percolation problem defined by 
the blockages is short-range correlated: given that a site 
is a blockage, its neighbors are slightly more likely to also 
be blockages. This detail shifts the percolation thresh¬ 
olds in d > 1 but does not modify the universal critical 
properties. 

The detailed thermalization properties of the phases 
in the semi-classical limit are revealed by the growth of 


FIG. 4: Standard deviation of saturated entanglement en¬ 
tropy for perturbed Glifford circuit in d = 1 at p^ = 1. The 
crossover in finite-size flow at ^ ~ 0.37r indicates the transition 
into the delocalized phase at large 0. 

Transition away from the Clifford limit — In this sec¬ 
tion, we relax the restriction to the Clifford group in order 
to explore the delocalization transition with fq 0. This 
severely limits the analytic and numeric tools available 
for studying the dynamics: the Gottesman-Knill theorem 
no longer provides a route to efficient classical simulation 
as even the addition of an infinitesimal single qubit rota¬ 
tion to the set of Clifford gates makes the system univer¬ 
sal for quantum computation m- Thus, the arguments 
to follow are somewhat more heurisic. 

We begin by assuming that the Clifford localized phase 
persists in the presence of weak generic perturbations. 
We will numerically check this assertion for d = 1 below; 
there are also several recent works arguing for the exis¬ 
tence of Floquet MBL phases in other models [32h34] . In 
the perturbed case, the lOMs develop exponential tails 
with typical length fq > 0. Near the percolation transi¬ 
tion, an lOM associated with an ergodic puddle is of large 
underlying size and the quantum tails simply fuzz the 
boundary if fq see Fig. The dressed lOMs now 
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interact with one another according to where dg 

is the distance of closest approach of the puddles. If the 
typical interaction scale is smaller than the level spacing 
on each puddle 2“^^), where ^ is the num¬ 
ber of sites in a typical puddle, then the states on the 
two puddles cannot hybridize and the lOMs remain lo¬ 
cal. This suggests the MBL phase persists up to a critical 

given by: 

^ ^{Pb- as Pb 4 , pI, (1) 

where we have assumed that de ^ 0(1), as there is a fi¬ 
nite density of points on the lattice where the puddles 
nearly connect. The generic phase boundary has the 
shape shown in Fig. 

There are a few possibilities for the nature of the phase 
boundary at non-zero ^q. The boundary could be first 
order. Assuming that it is continuous, however, the line 
of critical points should lie in the generic Floquet uni¬ 
versality class with no symmetries. By analogy with the 
semi-classical limit of the integer quantum Hall transition 
mnig, we do not expect this to coincide with classical 
percolation, as at = 0. Indeed, quantum tunneling 
across saddle points is well-known to increase the corre¬ 
lation length exponent of the Hall transition [43]. We 
therefore conjecture that the same inequality holds for 
the correlation length exponents of the MBL transition. 
We leave the construction of a detailed generic critical 
theory for future work. 

It is clear that the delocalized non-ergodic phase re¬ 
mains delocalized under generic perturbations. However, 
it is likely to become fully ergodic. Any local lOM in¬ 
teracts with the infinite percolating ergodic cluster with 
some finite strength for any > 0. However, the level 
spacing on the infinite cluster vanishes. Thus, by the 
arguments above, every local lOM hybridizes with the 
percolating cluster and delocalizes. 

Finally, we provide evidence in d = 1 that (i) localiza¬ 
tion persists away from the Clifford limit, and (ii) there is 
a transition at > 0. Numerically exact time evolution 
is limited to small system sizes N = 6 — 12; we therefore 
work deep in the localized regime at = 1. A layer of 
random single qubit rotations provides the perturbation 
to the Clifford brick-layer circuit. The rotation on qubit 
i is given by Vi = where a is the vector of 

Pauli operators, n is a unit vector on the Bloch sphere 
with polar angle 0 and azimuthal angle 7r/4, and 


is uniformly chosen. 

This particular model for quenched disorder can be 
motivated as follows. Recall that [Zi^U] = 0 for all i at 
P 5 = 1 in the Clifford limit. The single-qubit rotations 
dress the lOM Zi. If the axis of rotation is slightly tilted 
from the ^ axis by an angle d, we expect the dressing to 
be quasi-local; this leads to the generic MBL phase. At 


larger 0 > Oc, the dressing is non-local and the circuit 
has no local lOM. This should be the delocalized ergodic 
phase. 

As before, we perform global quenches from the anti¬ 
ferromagnetic state in the 2 :-basis and study the disorder- 
averaged saturated entanglement entropy [5 'a(oo)] be¬ 
tween the two halves of the chain after t ^ 2 ^^ time steps. 
We find that [5'y^(oo)] is much smaller than Na = N/2 
at small 0 while it approaches N/2 — c near 0 = nj2 (not 
plotted). This latter scaling is consistent with a fully er¬ 
godic volume-law; the constant correction c ^ log 2 e as 
proven in |44|. However, [5'y^(oc)] increases with N in 
both phases and the final state entropy density is not a 
sensitive enough measure at these system sizes to differ¬ 
entiate them. 

The standard deviation of 5 'a(oo), 

A5 = y[(<55^M)2], (3) 

however, behaves quite differently in the localized and 
delocalized phases m- In the MBL phase, the final en¬ 
tanglement entropy after the quench comes from the de¬ 
phasing of each of the Na local lOMs within region A. 
As each contributes an 0(1) number to Sa{oo), central 
limiting behavior implies that AS ^ \/N. In the delocal¬ 
ized phase, however, eigenstate thermalization suggests 
that there should be only exponentially small fluctua¬ 
tions of Sa with system size [45l [46]. Figure. shows 
the finite-size flow of this quantity as a function of 
and is clearly consistent with a transition at Oc ~ O.Stt, 
although the limited sizes prohibit systematic extraction 
of the N dependence of the flow. 

We note that a localization length ^q can be extracted 
from the residual polarization in the same experiment 
(not plotted), but its analysis is plagued by the same 
finite-size effects as the study in Ref. [23] . At the smallest 
value of O/t: = 0.05, the finite size flow is not visible 
within error bars because [£'^( 00 )] is itself very small. 

Conclusion — Periodic quantum circuits provide a use¬ 
ful perspective on localization. They distill the un¬ 
derlying entanglement transition from the complicated 
epiphenomena associated with globally conserved cur¬ 
rents [T31I17HS1]. Here, we identify disordered Clifford 
circuits as a tractable semi-classical limit with a delocal¬ 
ization transition that maps precisely to percolation, a 
well-studied beast [53]. The percolation exponents show 
up in dynamical properties such as the growth and satu¬ 
ration of entanglement entropy after a global quench. We 
conjecture they also lower bound the correlation length 
exponent in the generic case with quantum tunneling, 
just as in the integer quantum Hall plateau transition. 
To prove such a result, it would be very interesting to 
develop a network model for the interacting transition, 
perhaps in analogy with the recent work in one-dimension 
na. It would also be quite interesting to reintroduce 
global conservation laws, as these would likely lead to 
distinct universality classes for the delocalization transi¬ 
tion. 
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